Phase Transitions In Two Planar Lattice Models And 



Topological Defects: A Monte Carlo Study 

Subhrajit Dutta and Soumen Kumar Roy 
February 2, 2008 

Department of Physics, Jadavpur University, Calcutta— 700 032, India 

Abstract 

Monte Carlo simulation has been performed in the planar P2 and P4 mod- 
els to investigate the effects of the suppression of topological defects on the 
phase transition exhibited by these models. Suppression of the 1/2-defects 
on the square plaquettes in the P2 model leads to complete elimination of the 
phase transition observed in this model. However in the P4 model, on sup- 
pressing the single 1/2-defects on square plaquettes, the otherwise first order 
phase transition changes to a second order one which occurs at a higher tem- 
perature and this is due to presence of large number of 1/2 -pair defects which 
are left within the square plaquettes. When we suppressed these charges too, 
complete elimination of phase transition was observed. 

1 INTRODUCTION 

It is well known that conventional long range order cannot exist in a two dimensional 
continuous spin svstemlfTll. However the existence of topological charges leads to a 
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Table 1: The parameters obtained in the L=80, P4 model for the best fit g(r)= ar"^ + b of the 
correlation function at the temperatures indicated. In the normal lattice at these temperatures 
g(r) decays exponentially to zero. 

quasi long range order (QLRO) to disorder phase transition. The most notable and thor- 
oughly investigated example is the two dimensional XY model (0(2) model) in which, 
using a renormalization group technique, Kosterlitz and Thouless Q predicted a QLRO- 
disorder phase transition which is associated with the unbinding of the vortex-antivortex 
pairs (topological charges of strength ±1) which are stable topological defects in this 
system. The phase with QLRO is characterized by an algebraic decay of the spin-spin 
correlation function which is a slower decay than the fast exponential one which is ob- 
served in a completely disordered system. On the other hand in the two dimensional 
nonlinear sigma model, an example of which is the planar 0(3) model, there exists no 
stable topological defect and the system remains disordered at all finite temperature 0. 

Another class of the two dimensional systems of interest is the planar Lebwohl 
Lasher (LL) [4| model and a modified version of it to be elaborated below. In the 3- 
dimensional version of the LL model, the spins (of dimensionality 3), located at the 
sites of a simple cubic lattice, interact with nearest neighbors via a potential —P2{cos9) 
where 9 is the angle between the spins and P2 is the second Legendre polynomial. This 
model successfully describes the orientational aspects of a nematic and undergoes a 
weakly first order phase transition, representative of the nematic— isotropic transition, 
seen in a real nematic. A number of investigators 01611 have used a modified version 
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of the LL model by adding a — P4 temi to the usual — P2 one, P4 being the fourth or- 
der Legendre polynomial. The introduction of the P4 term reduces the sharpness of the 
peak of the P2 term in the potential at 6* = 7r/2 and may lead to the appearance of a local 
minimum, depending on the relative strengths of the P2 and P4 terms in the potential. 
This is found to enhance the first orderdness of the N-I transition. 

The two dimensional version of the LL model and a modified version of it with 
a pure P4 interaction between nearest neighbor spins have recently been investigated 
using Monte Carlo(MC) methods by a number of authors |I71[8l|9|. In the rest of this 
paper we shall refer to these as the planar P2 and P4 models respectively. Both models 
possess in addition to the usual 0(3) symmetry, a Z2 symmetry as well and this leads to 
the identification of the antipodal points in the order parameter space S^. The planar P2 
model is known to exhibit a continuous phase transition at a dimensionless temperature 
whose thermodynamic limit is 0.547 jSjl and the P4 model is characterized by a strongly 
first order transition at temperature 0.376 |9|. In the low temperature ordered phase in 
both models the pair correlation function shows an algebraic decay to a plateau which 
changes over to an exponential decay in the neighborhood of the phase transition|l8l|9|. 

The role of defects in the phase transition of various three dimensional spin systems 
is very difficult to study theoretically due to the nonlinearity introduced in the three 
dimensional nature of the spins. However simulation technique may be used to investi- 
gate the role of topological defects in these systems. Lau and Dasgupta lfTOl have shown 
numerically using the conventional Metropolis algorithm that monopoles (hedgehogs) 
are necessary for the phase transition in the three dimensional Heisenberg model. These 
authors observed that if one suppresses the formation of these defects in the 3-d Heisen- 
berg model, the system remains ordered at all temperatures and the transition to the 
disordered phase disappears altogether. The present work, which involves an elaborate 
MC study, was undertaken to investigate the effect of the suppression of the topolog- 
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ical defects on the phase transitions which the planar P2 and P4 models exhibit. The 
work was motivated to a large extent by the work of Lau and Dasgupta lfTDIl in the 3-D 
Heisenberg model. Another work, along the same line, which must be mentioned in this 
context, is that of Lammert et. al.| 1 1 1 who, in a MC study, have shown that the nature of 
the nematic-isotropic transition in a 3D nematic changes when one suppresses the for- 
mation of the stable line defects, called the disclination lines. Our work shows that the 
topological defects play a very crucial role in the phase transition in the planar P2 and 
P4 models and although these models possess the same symmetry they have remarkably 
different critical behavior. 

In the next section we briefly discuss the nature of the topological defects in the 
planar P2 and P4 models and the algorithms for their identification are presented in the 
following one. The details of our MC simulation is then presented followed by the result 
and discussion. 

2 Topological defects in the planar models 

In the two dimensional Heisenberg model there exists no stable topological defects. So 
the first or the fundamental homotopy group is just the set containing identity. However 
in the present planar models (P2 and P4 ) due to the local Z2 symmetry in addition to the 
0(3) symmetry there arises stable topological point defects known as 1/2-disclination 
points, where the director rotates through an angle of 180° around the defect core. The 
order parameter space is just the unit sphere with antipodal points identified. Any 
mapping of the other half integral point defects on the order parameter space is homo- 
topically equivalent to the mapping of the 1/2-defect. Point defects of integral strength 
are not stable in these models because of the so called 'escape to the third dimension'. 
Any attempt to escape from a configuration containing a 1/2 point defect leads to a more 
singular semi infinite line defect extending from the defect core. So the fundamental or 
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the first homotopy group of the concerned models is just the two element group Zo lfT^ . 

7ri = {0,l} 

It is known that topological instability does not necessarily imply physical instability O. 
If the path connecting the singular to nonsingular configurations of the free energy in- 
volves a configuration of higher free energy than either, then one may say that the topo- 
logically unstable singularity may possess a considerable degree of physical metasta- 
bility. This seems to happen in the P4-model and may be briefly explained as follows. 
Consider a configuration where each of the four spins at the lattice sites which form a 
unit square are in a plane and oriented at right angle to their neighbors. The —Pi{cos9) 
potential, besides having the global minimum at =0 (or 6=tx), also has a local mini- 
mum at = 7r/2. If the orientation of the spins are now gradually changed in order to 
make ^ (so as to reach the ground state ) a potential barrier will have to be over- 
come and the process becomes energetically costly. Thus there may exists metastable 
integral point defects in this model ifT^ . 

The algorithm for the detection of these defects is however nontrivial as it is not re- 
ally possible to enclose a 1 -defect by four spins alone. The possible method of detection 
of these defects is discussed in the next section. We add that we were unable to detect 
any such defect because of the low probability of their formation. 

3 The Defect finding algorithms 

In order to detect the 1/2 -point defects we have followed the algorithm originally pro- 
posed by Vachaspati lfT4l and subsequently used by others[ 15|. In order to trace out the 
topological defects it will be useful to see when a closed loop in the physical space will 
enclose a 1/2 disclination point. Let us consider a triangular plaquette ABC in the phys- 
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ical space. Due to local inversion symmetry we have to assign antipodal pair points on 
the unit sphere (S^) for each site of the triangular plaquette. Let points corresponding to 
A, B and C be (N,S) , (P,P^) and (R,R^) [ Figl]. Let us start from the north pole N. The 
point P or P^ is selected depending on which is closer to N. Let P be the selected point. 
So the arc NP of the great circle on the order parameter space is traversed when we go 
from site A to B on the physical space. Select from R or R^ whichever is closer to P. If 
the selected point is closer to S then the mapping of ABC is a non contractible loop and 
the plaquette will enclose a 1/2-disclination point defect. However if the selected point 
for C is closer to N, then the mapping is contractible to a point on the sphere and no 
defect will be enclosed by the triangular plaquette. It may be noted that if Si, Sj and S^ 
are the spin variables associated with the points A, B and C then the triangular plaquette 
will enclose a 1/2-defect if 



Sgn[{Si, Sj){Sj, Sk){Sk, Si)] = -1 (1) 

Priezev and Pelcovits Q in their work on 3-dimensional nematics have defined de- 
fect counting operators base on this principle. Beside these two mathematically equiva- 
lent methods, an algorithm for detecting 1/2-defects in RP^ models can be developed on 
the method first proposed by Berg and Luscher llT6ll . The method works as follows. The 
projection matrix P associated with each unit spin vector S in RP^ model obeys relation 
P^ =P and Tr P=l and its elements may be defined as Pa/3=S'aS'/3 where a ,(3 =1,2,3- The 
charge at a lattice site x* enclosed by an elementary triangular plaquette which has the 
projection matrices Pi , P2 and P3 associated with its corners is given by 



We have used and checked that all the three above mentioned algorithms for the 
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detection of the 1/2-defects in triangular plaquettes are exactly equivalent in all cases in 
both models. 

The detection of the metastable 1 -defects in the P4 model is a non-trivial job and the 
probability of their formation is very low as this requires a very 'special arrangement of 
the order parameter over many uncorrelated domains' 1 17 1 . In 3-d nematics, in principle 
where both line and point defects (hedgehogs) may form, no point defects are observed 
in experiments on quenched nematics ifTSi until the defect network has coarsened appre- 
ciably. It has been observed that the monopoles were only formed by string interactions 
and none were generated during the quench. Using the topology, more specifically, 
the homology of the order parameter space lfTTI . Hindmarsh has explained why the ex- 
pected density of point defects in extremely low. The observation, briefly speaking, is 
that in order to cover RP^ space twice (which is necessary for a topological 1 -charge) a 
roughly spherical arrangement of a minimum of twelve uncorrelated adjacent domains 
is necessary and this has a probability -10^®. In a two dimensional nematic the same 
consideration are believed to apply for the the -i-l point defects lfT9ll . We have used the 
algorithm using 12-spin configuration proposed by Zapotacky et. al. fT5\ to detect the 
1 -charge but could find none. 

4 The Simulation Details 

In the present paper we have used the conventional Metropolis spin update algorithm ll20l 
ED with periodic boundary condition in order to study the role of the topological de- 
fects in the phase transitions exhibited by the planar P2 and P4 model. Lattice sizes 
ranging from 20x20 to 80x80 were used and a part of the work was performed using the 
histogram reweighting technique of Ferrenberg and Swendsen ll22t . 

In order to carry out the procedure of the suppression of the defects in the planar lat- 
tice models we have included a chemical potential term associated with the topological 
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Figure 1: The Triangular plaquette ABC in physical space (left) and the order parameter space 
(right). On (N,S) is assigned for A, (P,P^) for B and (R,R^) for C. The great circle xy is 
perpendicular to NS and xiyi is perpendicular to PP^. If the point R is outside the region 
enclosed by two great circles then the corresponding loop is contractible, otherwise ABC will 
enclose a 1/2-disclination point, [see ref.l 14,11 

charges ifTOll . The Hamiltonian used in the simulation is given by, 

H = -J2Pl {cosOij) + A ^ Qijki (3) 

ij ijkl 

where L is either 2 or 4. Oij is the angle between the nearest neighbor spins i, j and 
the Qijki is the sum of charges of two triangular portions of a square plaquette. A 
positive A makes the formation of the charges expensive in terms of energy and for 
almost total suppression of the charge a large value of A (about 10 to 60, but independent 
of temperature) was normally chosen. In order to obtain the unrestricted simulation we 
set A = 0. The charge enclosed by ijk for instance is given by. 

Qijk = ^[1 - sgn{{Si, Sj){Sj, Sk){Sk, Si)}] (4) 

Clearly the sum Qijki can be 0, 1/2 or 1. If Qijki is then the square plaquette encloses 
no charge. If it is 1/2 then a 1/2-disclination point is enclosed. But if it is 1, then it 
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should not be confused with an integral point defect. In fact this corresponds to two 
closest possible 1/2-charges situated within a square plaquette of linear dimension 
equal to lattice spacing. 

In both the planar models we have investigated, the smallest part of the system in 
real space that can enclose a 1/2-point defect is a triangle. Each elementary square 
plaquette can be diagonally cut into two triangles and if these two adjacent triangles 
each enclose a 1/2-defect then this leads to a 1/2-pair charge being enclosed by the 
elementary square. If only one of those triangles encloses a 1/2-charge then the square 
in turn encloses a single 1/2-charge. We denote number of square plaquettes enclosing a 
pair of 1/2-charges by ni (this should however not be confused with a topological defect 
of charge 1). Similarly the number of elementary squares enclosing a single 1/2-charge 
is denoted by ni/2. In a recent work, Mondal and RoyfTSl have observed that the ratio 
ni/2/ni behaves like a response function in both the models although their behavior is 
different in the two systems. When plotted against temperature, the ratio ni/2/ni exhibits 
a maximum at the transition in the P2 model while in the P4 model it shows a sharp fall 
at the transition. Finite size effects are also prominent in the transition temperatures 
thus obtained from the ni/2/ni vs T plots. 

In our MC simulation, while investigating the effect of suppression for the charges in 
the two planar models, we have treated the single 1/2-charge and the 1/2 -pairs (within an 
elementary square plaquette) on different footings. For the simulation where no charge 
suppression has been attempted we have set A = 0. For the suppression of the charges 
represented by ni/2 , A 7^ only for Qijki =1/2 while for total (both single- 1/2 and 
1/2-pair suppression), A 7^ for Qijki 7^ 0. In the P2 model, complete suppression 
of single 1/2-defects was found to lead to complete suppression of the 1/2-pair defects 
and this leads to complete elimination of the phase transition in this model. In the P4 
model, however, the suppression of the single 1/2-defects leaves a large number of 1/2- 
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pair defects within the elementary squares and the evidence of a new phase transition 
at a higher temperature is obtained. When these defects too were suppressed, the phase 
transition totally disappears. 

In order to estimate the critical exponents and the thermodynamic limit of the critical 
temperature of the new transition which we obtained in the P4 model and which seems 
to be of second order, we have applied the finite size scaling method. Finite size scal- 
ing method is a technique of estimating the critical exponents and the thermodynamic 
limit of the transition temperature by observing how the measured quantities vary with 
the system size. In finite size scaling method (the data collapse method in particular) 
we extract the part of the thermodynamic function which does not contain system size 
explicitly |l2T|. This part is called the scaling function. If proper values of the critical 
exponents and the thermodynamic limit for the transition temperature are chosen then 
the scaling function for different system sizes get collapsed. In this paper we have used 
the data collapse technique for estimating the critical exponents associated with the spe- 
cific heat and the order parameter. The critical exponents associated with the correlation 
length, specific heat and the order parameter are denoted by u, a and (3 respectively. 

For specific heat, the scaling relation stands as 

Cv = L^^'^C^L^t) (5) 

where t is the reduced temperature and C is the specific heat scaling function. 
Similarly for the order parameter ( <q^>)(see eq. |71), the scaling relation is given by 

< >= L-^^/'QiLh) (6) 
where Q is known as the order parameter scaling function ifTDl . 
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5 RESULTS AND DISCUSSION 



We have evaluated various thermodynamic properties like internal energy per particle 
(<E>), specific heat, order parameter etc. The specific heat was evaluated by taking the 
temperature derivative of the <E> as well as from fluctuation of the energy. Due to the 
local inversion symmetry, the order parameter is a second rank tensor. 

As a measure of the order prevailing in the system we used the quantity <q^> given 
inref.irm 

where Qab = ■^J2Qi'>')ab^s the nematic tensor order parameter, where 
Q{'i)ab = {naUi, — ^5ab) {n is the molecular axis of the i*'^ molecule) and N is the total 
number of sites (i) in the lattice. This definition ensures that <q^> is zero in a fully 
disordered system and 1 for a fully ordered system. 

In case of P2 model we have simulated for linear dimension L=40 and 60. In Fig2 
we have depicted the specific heat versus temperature plot for the P2 lattice model of 
size 40x40. The unrestricted simulation shows a peak which disappears when the single 
1/2-charges on the square plaquettes were suppressed. We have observed that the after 
suppression of the single 1/2-charges on the square plaquettes there remains no 1/2 -pair 
defect in any square plaquette. The temperature dependence of the order parameter used 
in this model for the two cases are shown in Fig3. While the temperature derivative of 
<q^> has a peak at the transition temperature in the normal case (where no defect is 
suppressed), which presumably is a signal of a phase transition in a finite system ll23ll . 
in the defect free case it seems to lose the characteristic shape and shows a smooth and 
rather slow decrease with temperature and vanishes at around a temperature T=6. We 
would be inclined to conclude from the results on Cy and <q^> that the defect free phase 
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exhibits no phase transition at all. 

Turning to the P4 model, which has a characteristic strongly first order phase transition^, 
we first point out that the suppression of the single 1/2 defects on the square plaquettes 
here does not result in suppression of the 1/2-pair defects on the square plaquettes. This 
observations is different from what happens in the P2 system where the suppression 
of the single 1/2 defects on the square plaquettes leads to suppression of the 1/2-pair 
defects. On suppressing the single 1/2-defects a new phase transition is observed. We 
however point out that it is impossible to make the system completely free of topological 
defects, even when arbitrarily large values of A are used. However the residual charges 
left were of very insignificant amount. For instance, at T=0.55, the traces of the sin- 
gle 1/2 and 1/2-pair charges which could not be suppressed were about 0.02% of these 
charges present in the system at the same temperature after single 1/2-charge suppres- 
sion. Fig4 shows the temperature dependence of the specific heat for the normal lattice 
(L=40) and after suppression of the single 1/2-defects ( for L=40 and 60). The peak in 
Cy for the latter case, while greatly reduced in size (from 80 to about 17), shifts to a 
higher temperature which for both lattices is close to 0.494. Presumably the phase tran- 
sition which we observed after suppressing the single 1/2 defects is due to the presence 
of the 1/2-pair defects. When we suppressed the 1/2-pair defects too in all triangular 
plaquettes of the lattice no evidence of any peak in Cy was observed, although a hump 
like feature was seen with Cy -3.5 over a temperature range extending from 0.53 to 
0.62. In case of P4 model we have also investigated the response of the defect density to 
the phase transition. When no kind of charge suppression was applied a large number 
of single 1/2 defects as well as 1/2-pair defects were found to be present. In fig5 the 
temperature dependence of the density of single 1/2-defects are shown for unrestricted 
simulation. In fig6 the temperature dependence of the density 1/2-pair defects is shown 
for both before and after suppression of the single 1/2 defects. Clearly the temperature 
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Figure 2: The specific heat, vs. temperature plot in P2 model. + corresponds to Cy obtained from 
energy fluctuation in the L=40 lattice and the continuous curve is taken from refjH where it has 
been obtained using multiple histogram reweighting and peak is at T=0.587. Both corresponds 
to normal MC simulation with no charge suppression. The x represents Cy for L=40 lattice after 
suppressing the single 1/2 defects on the square plaquettes. 

derivative of defect density behaves like other response functions for all the cases. 

In figure 7, we have depicted the temperature dependence of the order parameter 
<q^> for all three cases in the P4 model for L=40. The peaks of temperature derivative 
of <q^> occur respectively at 0.379 (normal MC) and 0.490 (after suppressing single 
1/2-defects) and their heights are 1 19.5 and 18.5 respectively. On suppressing all the 1/2 
charges on all triangular plaquettes of the lattice (total charge suppression), d<q^>/dT 
no longer exhibits a peak and <q^> seems to have an asymptotic value of 0. 1 at T=2, 
up to which the investigation has been made. We add that d<q^>/dT after total charge 
suppression has a feature similar to that seen in Cy in that a broad peak of height ~4 was 
seen over a temperature range from 0.52 to 0.6. The existence of this broad hump of 
insignificant magnitude over an extended temperature range can not be a sign of phase 
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Figure 3: The order parameter <q^> obtained for the L=40 lattice in the P2 model with no charge 
suppression. The temperature derivative has a peak at T=0.567. The inset shows the <q^> vs T 
plot after suppressing the single 1/2-defects on the square plaquettes. 
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Figure 4: Cy plotted against T for unrestricted simulation of L=40, P4 lattice (left) and the 

same for L=40 and L=60 lattice obtained after suppressing the single 1/2-defects on the square 
plaquettes. The slightly bigger peak is for the L=60 lattice. The peak occurs at T=0.494. 
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Figure 5: The density of the single 1/2-defects enclosed by the square plaquettes in case of L=40 
P4 model with no charge suppression. 
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Figure 6: The density <ni > of 1/2-pair charges for unrestricted case of L=40 P4 model (left) 
and the same after suppressing the single 1/2-defects . The temperature derivative of the former 
has a peak at T=0.379 and the latter(inset) peaks at T=0.494. 
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Figure 7: The temperature dependence of <q^ > for the L=40 P4 model. On the left is the 
unrestricted MC result and on the right is the case when the single 1/2-charge on square plaquette 
were suppressed. The inset shows <q^> when all charges are suppressed. The temperature 
derivative of the first two curves have peaks at 0.379 and 0.490 respectively. 
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Figure 8: The correlation function g(r) plotted against r for L=40 P2 lattice at T=0.6 for both 
unrestricted case and the suppressed case. 

transition. These, in someway, may be connected to the existence of the small number of 
residual charges left in the system after the attempt to suppress them completely failed. 

We now turn to the pair correlation function g(r)=< cos^9ij—l/?> > where r is 
the separation between the two spins i and j which make an angle 6ij with each other. 
This function for the P2 model at T=0.6, for instance, which is a temperature much 
higher than the normal critical temperature in this model, decays exponentially to zero, 
as it should, in the complete absence of long range order and a best fit with g(r)= 
aexp{—\r) yields a=0.439 and A=0.248. For the single 1/2-charge (and consequently 
of 1/2-pair charges) suppressed case, it decays algebraically and a best fit like g(r)= 
ar-P + 6, yields the parameters a=0.404, b=0.022 and p=0.372 [Fig 8]. In figure 9, the 
g(r) vs r plots for the P4 model are shown for T=0.48 (which is greater than the normal 
transition temperature but less than the transition temperature obtained after 1/2-charge 
suppression), T=0.5 (which is slightly higher than the observed transition temperature 
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Figure 9: The correlation function g(r) plotted against r for the L=80, P4 lattice. For the curves 

(a), (b), (c), T=0.68, for the curves (d), (e), (f), T=0.5 and for (g), (h), (i), T=0.48. The three 
curves in each set correspond to the normal MC, 1/2-charge suppressed and total charge sup- 
pressed cases. The parameters used to fit these curves are listed in Table 1. 
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Figure 10: Free energy vs energy for P4 model for L=40(upper) and L=60(Iower) sizes after 
suppressing the single 1/2-defects on the square plaquettes. 



after 1/2-charge suppression) and T=0.68. At all of these temperatures, as one would 
expect, the correlation function in the normal case decays exponentially to zero. The 
two other cases give best fits for an algebraic decay to a plateau, and the parameters 
are listed in the Tablel. We find that in the P4 model the results of single 1/2-charge 
suppression is the same as that in the P2 model while with the suppression of the 1/2-pair 
charges too, the asymptotic value of the order prevailing in the system increases further. 
The phase transition, which the P4 system has after suppressing the single 1/2-charges 
is second order. We have evidence that it is not first order as the dual peak structure 
of the probability distribution as a function of energy, which is so distinct in the P4 
model[9|, has been found to disappear totally after suppressing the single 1/2-charges 
on the square plaquettes. In fig 10 we have given the the free energy of 40x40 and 60x60 
lattice in P4 model after suppressing the single 1/2 defects on the square plaquettes. The 
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Figure 11: The order parameter <q^> plotted against temperature T for three lattice sizes indi- 
cated for the P4 model after suppressing single 1/2-defects. 
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Figure 12: The specific heat Cy vs T for the lattice sizes indicated for the P4 model after the 
single 1/2-charge suppression. The values of Cy here obtained from d<E>/dT and the histogram 
reweighting technique l22ll was used. 
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Figure 13: Collapse of the order parameter for L=20,40,60 and 80 sizes. The best collapse 
obtained at Tc(oc) = 0.498 (thermodynamic limit), v =0.99 (exponent for correlation length), 
/? =0.16 (exponent for order parameter), t is the reduced temperature, (T-Tc) /Tc 
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Figure 14: The collapse of the specific heat data for system sizes L=40, 60 and 80. The best 
collapse obtained at thermodynamic limit of the transition temperature (Tc(oo)=0.4938), v = 
0.94 and a(exponent for the specific heat) = 0.077. 
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Figure 15: Evolution of order parameter for 80x80 P2 model after suppressing the 1/2-defects 
enclosed by the square plaquettes using A =60. The final values of order parameter for three 
initial states with qo^ =1.0, qo^ =0.19 and qo^=0.09 are almost same. 
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Figure 16: Evolution of total energy of the 80x80 P4 model after suppressing the 1/2-defects 
enclosed by the square plaquettes using A=60. The final values of Energy for two different 
initial states are almost same. 
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single well structure of the free energy indicates the second orderdness of the new phase 
transition after suppressing the single 1/2-defects. 

We have also used standard finite size scaling method available for second order 
phase transition lOn in order to estimate the critical exponents and the thermodynamic 
limit of transition temperature of the new phase transition that was observed after sup- 
pressing the single 1/2-defects. In Figll and Figl2 we have depicted the order pa- 
rameter <q^> and the specific heat Cy plots respectively as function of temperature for 
different lattice sizes after the single 1/2-charge suppression in the P4 model. We have 
used standard data collapse technique to collapse the data of the figures 11 and 12 and 
the resulting diagrams are Fig 13 and Fig 14 respectively. In fig 13 we have shown the 
collapse of the order parameter scaling function for L=20, 40, 60 and 80. While in fig 14 
the collapse of specific heat scaling function for the system sizes 40,60 and 80 are dis- 
played. The data collapse clearly shows the phase transition that P4 model exhibits after 
the single 1/2-charge suppression is second order. The parameters we have obtained are 
Tc=0.498, (3 =0.16 and u=0.99 for the collapse of order parameter and =0.4938 , 
a =0.077 and u =0.94 from collapse of the specific heat. We would remark that the 
high temperature behavior of the order parameter <q^> after total charge suppression 
does not indicate the presence of a phase transition. This is true for both models. 

We have used large values of A in our simulation in order to suppress the evolution 
of defects in both models. It is known that any Monte Carlo study is faithful only if we 
can reach any point in the phase space starting from any other point . So there must be 
a path connecting the two points in phase space with non zero probability. This actually 
indicates that we should be careful that we are not trapped in any small region of the 
phase space. We have investigated the phase space connectivity in both the models by 
observing the evolution of order parameter or energy with MC steps. The connectedness 
is satisfied if the observed quantities for different initial states converge to the same final 
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value. In fig 15 we have shown that in case of the 80x80 P2 model, after suppressing the 
1/2-defects ( by using A=60 ) on the square plaquettes, the final values of the order 
parameter is same for three different initial configurations. Similarly in Figl6 we have 
shown that the same thing happens for the total energy of the 80x80 P4 model after 
suppressing the 1/2-defects on the square plaquettes. It is therefore clear that we can 
use a value of A at least up to 60 without violating the phase space connectivity. 

6 CONCLUSION 

It is established in this paper that topological defects play a very important role in the 
phase transitions exhibited by the two planar lattice models we discussed. The observed 
difference in the critical behavior seems to be due to the difference in the role played 
by the topological defects. It is shown in this paper that for the phase transition in both 
the models topological defects are necessary. In the P2 model the phase transition is 
governed by the 1/2-disclination points enclosed by the square plaquettes only (single 
1/2). On suppressing these single 1/2-defects we have shown that the phase transition 
was totally eliminated. However the picture is different in case of two dimensional 
P4 model. Where on suppressing the single 1/2-defects the nature and the transition 
temperature of the phase transition gets changed. We have also shown that the new phase 
transition in the P4-model is due to the presence of large number of 1/2-pairs (enclosed 
by two triangular portions ) within the square plaquettes which remains unsuppressed 
even after suppressing the single 1/2-defects. This leads to another important conclusion 
that for tracing out all the disclination points in two dimension, the minimum closed 
loops in physical space to be considered are the smallest triangular cells formed by 
three nearest neighbor sites. 

It will be interesting to apply the idea of triangular plaquettes in three dimensional 
nematics, where the stable topological defects are disclination lines as well as monopoles. 
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Let us consider small disclination loop in a three dimensional nematic [fig 17] . In the fig- 
ure a cubic unit cell ABCDEFGH is shown. The disclination loop has crossed the face 
ABCD twice. However if we take mapping of ABCD on the order parameter space then 
we would get a contractable loop in the order parameter space and conclude that ABCD 
does not enclosed any line defect. It is true that ABCD does not enclose a single line, but 
it encloses two close lines (as the loop intersects the face ABCD twice), which could be 
detected only if we consider the two triangular portions (ABC and ADC) of ABCD. So 
considering the elementary plaquettes to be square plaquettes of linear dimension equal 
to lattice spacing, we would miss loops as shown in the figure. In real nematics small 
loops are known to carry monopole (hedgehog) charges. Actual monopole structure is 
very rare in real nematics. As already stated, the low probability of actual monopole 
structure is discussed by Hindmarsh |17|. The monopoles that comes in to play in the 
three dimensional nematics arises mainly from the small disclination rings. In our work 
it is evident that in P4 model it is very important to consider these triangular plaquettes 
in order to trace out the topological defects properly. If the same work is carried out on 
three dimensional P4 model and disclination lines were traced out using triangular pla- 
quettes, then one is likely to get a large number of rings as shown in the figurelV. The 
small disclination rings which behave like monopoles, is expected to have important 
role in the phase transition of the three dimensional P4 model. 
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Figure 17: The cubic unit cell ABCDEFGH in a three dimensional nematic. The disclination 
ring cuts the upper face ABCD at X and XI. The triangular portions ABC and ADC must 
enclose a 1/2 -defect each. 
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